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^ Abstract. A high-resolution A^-body simulation suggests that stellar 

streams in the discs of barred galaxies are common and strongly time- 
CJ . dependent. The velocity distribution of stars in the Solar neighbourhood 

^^ ' betray many such streams, including a stream of outward moving stars 

with low angular momentum. This stream is interpreted as a signature 
of the Galactic bar, in the sense that its stars have just enough energy 
(Jacobi's integral) to cross the corotation resonance. 
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The Herculis stream 



0> ■ Figure || shows the observed velocity distribution in the U-V plane, based on the 

^^ \ Hipparcos input (radial velocities) and output (tangential velocities) catalogues. 

The distribution is characterised by several distinct streams, and in particular a 
Qh| stream with an asymmetric drift of 40 — 50 kms^^ and [/ > 0, which hereafter 

Q ■ we will refer to as the "Herculis" stream according to a comoving Eggen group 

H '. (Skuljan et al. 1999). This stream most probably has a dynamical origin because 

t/^ ■ its stars are mainly older than a few Gyr and present a large range of metallicities 

' (Raboud et al. 1998), and is under-represented in the Hipparcos sample, which 

is biased towards young stars. The average luminosity of stars in this sample 
►v> \ indeed increases with distance and the sample only covers the vertical region of 

\^ • the galactic plane where the fraction of young stars is largest. 
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2. Simulation 

We have realised a large 3D N-hody simulation designed to evolve into realis- 
tic barred models of the Milky Way. The initial conditions are adapted from 
Fux (1997) and N = 1.4 x 10^ particles, allowing to compute very detailed 
velocity distributions in small space regions. 

After the formation of the bar (at t ~ 1.4 Gyr), the simulation reveals a 
complex velocity structure with multiple streams occuring almost everywhere 
in the disc. The streams, as well as the velocity dispersions, remain very time- 
dependent until the end of the simulation (i = 2.65 Gyr), even within regions 
corotating with the bar. Outside corotation and for an inclination angle of the 
bar < (/3 < 90° , the simulation often displays a Herculis-like stream (Figure ||) . 
The model selected here is a snapshot at t = 2.15 Gyr, scaled to a corotation 
radius of {R12 + Ri5)/2 = 4.5 kpc, where Ru and R45 are the galactocentric 
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Figure 1. U — V velocity 
distribution of the Hippar- 
cos single stars with distance 
d < 100 pc, a(7r)/7r < 0.1 
and radial velocities in the 
input catalogue (3481 stars). 
The velocities are relative to 
the Sun; V is positive to- 
wards galactic rotation and 
U towards the anti-centre. 



distances of the Lagrangian points -^1/2 and L4/5, and to a circular velocity in 
the azimuthally averaged potential of 200 kms~^ at i?o = 8 kpc. In these units, 
RoLR = 7.7 kpc and ajj : ay '■ (Jw = 55 : 34 : 21 kms^"*^ at Ro and (p = 25°. 



3. Interpretation 

Our interpretation for the Herculis stream follows the same lines as in Raboud 
et al. (1998). In a rotating barred potential, the value of the Hamiltonian: 



H(x, x) 






(1) 



where $eff = ^{x) — ^^p{x'^ + y'^), $(x) is the gravitational potential and ilp the 
pattern speed of the bar, is an integral of motion (Jacobi's integral). Three crit- 
ical values of this integral are associated with stars corotating at the Lagrangian 
points {H12 = ^csiLi/2) and H45 = <I>eff(-^4/5)) and stars on circular orbits at 
the outer Lindblad resonance in the axisymmetriesed potential (-ffoLR)- Stars 
with H > H12 follow "hot" orbits which are susceptible to cross the corotation 
radius. In the simulation, even far beyond the solar circle, a substantial fraction 
of particles are on such orbits and the Hamiltonian distribution frequently shows 
a secondary peak between H12 and H45 (Figure ^). If U and V are measured 
relative to the galactic centre, Equation (||) also writes: 



{V -R- npf + U^ + W^ = 2{H - $cff). 



(2) 



Thus, neglecting the vertical velocity component W, the contours of constant 
Hamiltonian in the U — V plane are circles centred on {V, U) = {R ilp, 0) and of 
radius 2{H — ^cs)- Figure |2| shows that the Herculis-like stream in the model 
falls just outside the -fri2-contour and thus concerns particles on hot orbits. 
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Figure 2. Model U — V distribution as a function of position in the 
galactic disc (t = 2.15 Gyr). The corotation is at 4.5 kpc. The dashed- 
dotted, dashed and dotted curves are the contours of constant Hamil- 
tonian for H = H/^^, H12 and Hqlr respectively (assuming W = 0). 
The circles indicate the circular orbit in the axisymmetriesed potential. 
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Figure 3. Hamiltonian distribution of the disc particles at i? = 11 kpc 
{t = 2.15 Gyr). The dashed-dotted, dashed and dotted vertical lines 
respectively indicate the values H45, H12 and -f^oLR- The part of the 
distribution with H > H12 represents the hot particles which may enter 
inside corotation. 
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Figure 4. Left: model U — V distribution for R = 8 kpc and (p = 25° 
{t = 2.15 Gyr). The circular contours are as in Figure and the bold 
dots indicate the disc particles left of the diagonal line which stayed 
longer than 300 Myr inside corotation during the past Gyr. Right: 
orbits of these particles integrated 525 Myr backwards. The dotted 
lines are contours of constant ^^s and the pattern rotation is clockwise. 



This interpretation is different from the one proposed by Dehnen (1999), 
who attributes the valley between the Herculis stream and the Hyades-Pleiades 
branch to stars scattered off the OLR. In our case, the valley is produced by the 
decline of the hot orbit population as /7 — s- Hi2. But in both cases the position 
of the OLR is found near Ro, in agreement with Kalnajs (1991) conclusion based 
on the Hyades and Sirius streams. 

Figure Q traces back in time the orbits of the particles in the Herculis-like 
stream which were recently captured inside corotation (using the instantaneous 
frozen rotating potential). These particles come from the low Hamiltonian edge 
of the stream and represent about 10% of the particles in the stream. This 
percentage strongly depends on R and (p. 
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